Anomalous Paramagnetic Magnetization in Mixed State of CeCoIn 5 single crystals 
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Magnetization and torque measurements were performed on CeCoIns single crystals to study the 
mixed-state thermodynamics. These measurements allow the determination of both paramagnetic 
and vortex responses in the mixed-state magnetization. The paramagnetic magnetization is sup- 
pressed in the mixed state with the spin susceptibility increasing with increasing magnetic field. The 
dependence of spin susceptibility on magnetic field is due to the fact that heavy electrons contribute 
both to superconductivity and paramagnetism and a large Zeeman effect exists in this system. No 
anomaly in the vortex response was found within the investigated temperature and field range. 
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Introduction 



The interplay between superconductivity and mag- 
netism is a subject of great interest in the study of su- 
perconductors. The heavy fermion material CeCoIns is 
a strongly correlated /-electron superconductor, which 
makes it a good candidate to study this effect. This ma- 
terial displays several novel phenomena. For example, it 
is in the vicinity of the antiferromagnetic quantum crit- 
ical pointji^ As a result, its magnetic susceptibility x 
diverges at low temperature T as x k y-o.42 (R e f s> 4^ 
5). Heavy electrons are essential for the development of 
superconductivity**'^ Angular dependent thermal trans- 
port and specific heat measurements in a magnetic field 
provide evidence for d-wave pairing symmetry, which in- 
dicates singlet pairingi^iiS Nuclear magnetic resonance 
(NMR) measurements report suppressed spin suscepti- 
bility in the mixed state as a function of tempraturei* 
Large spin fluctuations exist in this system. There is an 
unusually large specific heat jump at the superconducting 
transition temperature T c q, which is due to the supercon- 
ducting pairing, but also to strong spin fluctuations! 11 ! 12 

A magnetic field suppresses superconductivity by cou- 
pling to either the spins or the orbits of the electrons. 
If the spin effect dominates, then the material is in the 
Pauli limit. On the contrary, if the orbital effect dom- 
inates, then the material is in the orbital limit. The 
Maki parameter a = y/2H C 2o/H p (where £f C 20 is the or- 
bital critical field in the absence of the Pauli limiting 
and Hp is the upper critical field limited by Pauli para- 
magnetism) gives the relative strength of the orbital pair 
breaking by magnetic field and Pauli limiting^ In the 
standard Bardeen-Cooper-Schrieffer (BCS) theory, the 
orbital effect dominates the Pauli limiting effect. This is 
the case for most superconducting materials. However, 
heavy fermion materials have large effective mass m*, so 
the Fermi velocity is very small, hence, the orbital effect 
is greatly reduced in heavy fermion materials. In partic- 
ular, CeCoIn 5 has a small Fermi energy, large supercon- 
ducting gap and a short coherence length. Also it is in 



the clean limit, with long mean free path, which is much 
larger than the superconducting coherence length. The 
value of Maki parameter a « 3.6. Hence, CeCoIns satis- 
fies all the theoretical requirements for the formation of 
the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) statej 14 i 15 
In addition, there is experimental evidence for the exis- 
tence of the FFLO state j 14 i 16 ' 17 All these findings support 
the fact that CeCoIns is, indeed, in the Pauli limit, which 
means that the spin effect dominates the orbital effect in 
this material at low temperatures. 

The usual orbital depairing effect forms vortices in 
CeCoIns in the presence of an applied magnetic field, 
while the Zeeman depairing effect forces the spins to 
align with the field, hence destroys the spin singlet pair- 
ing required for the existence of the Cooper pairs. For 
these reasons, one would expect an unusual mixed state 
for CeCoIn5, in which the diamagnetic and paramag- 
netic contributions could be anomalous in the presence 
of Zeeman effect. Therefore, it is important to address 
the issue regarding the paramagnetic and diamagnetic 
response in the mixed state and how does the Pauli 
paramagnetism affect the mixed-state thermodynamics 
of CeCoIns. This is essential to the understanding of the 
interaction between superconductivity and magnetism in 
heavy fermion materials. 

We performed magnetization and torque measure- 
ments in the normal and mixed states of CeCoIns in 
order to address the above issues. We successfully sep- 
arated the paramagnetic and vortex contributions. The 
paramagnetic magnetization is unusual and it has a non- 
linear magnetic field dependence, while the susceptibility 
X P in the mixed state increases with increasing field. The 
increase of the susceptibility with increasing field is due 
to the fact that heavy electrons contribute to both super- 
conductivity and paramagnetism and the Zeeman effect 
is large. The vortex contribution has no anomaly within 
the investigated temperature range, although Pauli lim- 
iting effect is present in this system. 
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B. Experimental Details 

High quality single crystals of CeCoIns were grown us- 
ing a flux method. To remove the excess indium left 
on the surface of the crystals during the growth process, 
the crystals were etched in concentrated HC1 for several 
hours and then rinsed throughly in ethanol. The mass 
of the single crystal for which data are shown here is 5.5 
mg and the zero-field superconducting transition temper- 
ature T c0 = 2.3 K. 

Both dc magnetization and angular dependent torque 
measurements were performed in normal and mixed 
states, over a temperature range 1.76 K < T < 20 K 
in magnetic fields up to 14 T. The dc magnetization 
measurements were carried out using a superconducting 
quantum interference device magnetometer in magnetic 
fields applied parallel to the c-axis of the single crystal. 
The torque measurements used a piezoresistive torque 
magnetometer. The single crystal was rotated in the ap- 
plied magnetic field between H || c-axis (9 — 0°) and 
H || a-axis (6 = 90°) and the torque was measured as a 
function of increasing and decreasing angle, under vari- 
ous temperature-field conditions. Details regarding the 
background subtraction in the torque measurements can 
be found elsewhere. 18 



C. Results and Discussion 

Figure 1 shows the field dependence of the measured 
magnetization M mes at 1.76 K for H \\ c— axis. This is 
a representative M mes (H) curve in the mixed state for 
temperatures up to 2.10 K. For each M mes (H) curve, we 
zero-field-cooled the single crystal to the desired temper- 
ature and measured the magnetization in increasing field 
up to 50 kOe and then decreasing the field to zero. Note 
that M mes (H) is irreversible in the low field region and 
it becomes reversible above a certain H value. Also, the 
magnetization increases monotonically with increasing H 
up to a certain value, beyond which it becomes linear in 
H. We define this latter H value as the upper critical 
field along the c-axis W^iT). 

Plotted in the inset to Fig. 1 is the H'^iT) phase 
boundary. The open squares are data taken from pre- 
vious reports^ while the open circles are data ex- 
tracted from the present measurements of M mes (H), 

with H^iT) defined as just described above. Note that 
all data fall on the same curve, which confirms that our 
definition of h]}£(T) is correct. 

Recently we reported large paramagnetism in the nor- 
mal state of this material^ As a result, the magnetiza- 
tion in the mixed state has two contributions: paramag- 
netic contribution and diamagnetic contribution due to 
the vortices; i.e., 

M mes = M p + M v . (1) 
Also, as discussed above, in the mixed state of CeCoIns 



one expects that orbital and Zeeman depairing mecha- 
nisms coexist. As a result, CeCoIn 5 could display a novel 
mixed state with anomalous paramagnetic and diamag- 
netic contributions. As a starting procedure to determine 
these contributions, we first assume the simplest case in 
which the paramagnetic magnetization is the same in the 
normal and mixed states, i.e., it is linear in H and the 
spin susceptibility is field independent. This has pre- 
viously been done in the study of Ndi.ssCeo.isCuC^-y 
(Ref. 20). Hence, we fit the linear part of the M mes {H) 
curve in the normal state and extrapolate it into the low- 
field region, where the sample is in the mixed state (see 
the solid line on the main panel of Fig. 1). By subtract- 
ing the paramagnetic magnetization in the normal state 
M n = as determined (Xn i s the normal state sus- 

ceptibility in the c direction), we should obtain the field 
dependence of the diamagnetic magnetization. 

The main panel of Fig. 2 shows the diamagnetic mag- 
netization Mi(H) = M mes (H) — x c n H & t 1-76 K and over 
the whole measured magnetic field range, i.e., from 
to 50 kOe, determined as just discussed above. Note 
that Mi(H) is non-monotonic with two peaks present 
in —Mi(H) curves: the first peak is at a very small field 
value (25 Oe for T — 1.76 K) and is very sharp. This peak 
corresponds to the lower critical field. A second, broader 
peak, however, appears at higher fields (for T = 1.76 K, 
this peak is in the field range 10 — 20 kOe). We show 
the enlarged non-monotonic part of the diamagnetic re- 
sponse Mi(H) in the lower inset to Fig. 2 for the mea- 
sured temperatures of 1.76, 1.80, 1.85, 1.90, 1.95, 2.00, 
2.05 and 2.10 K, from bottom to top. As the tempera- 
ture increases, the second peak becomes flatter, it shifts 
to lower H values, and at 2.00 K it disappears and the 
Mi(H) curve becomes monotonia Nevertheless, even at 
this temperature, the M\(H) curve does not resemble a 
typical diamagnetic M(H) curve. 

At a first glance, the second peak in —M\{H) looks like 
the second magnetization peak which appears in high- 
temperature superconductorS ) 22 i 23 or in the UPt3 heavy 
fermion material^ The reasons for the presence of the 
second magnetization peak in these materials are an en- 
hanced pinning and/or the presence of a phase transition. 
However, the second peak observed here in CeCoIn.5 is 
not due to enhanced pinning since M\(H) shows only a 
very small hysteresis even at the lowest measured tem- 
perature of 1.76 K. Also this second peak in —Mi(H) 
is not due to a phase transition since this peak is very 
broad. This indicates that the subtraction of the linear in 
H paramagnetic magnetization in the mixed state, which 
gives rise to the second peak in —M\(H) and which as- 
sumes that the paramagnetic magnetization in the mixed 
and normal states is the same, is not correct. So, we 
conclude that there is another contribution to the para- 
magnetic magnetization. Under these circumstances, in 
principle, it is very hard to separate the vortex and para- 
magnetic responses. However, we show here that torque 
measurements along with magnetization measurements 
permit the successful determination of both responses. 



3 



The above discussion, which points towards the pres- 
ence of another contribution to the paramagnetic mag- 
netization that is not linear in H, is consistent with the 
theoretical report of Adachi et al. which shows, based on 
quasi-classical Eilenberger formalism, that the functional 
form of the mixed-state paramagnetic magnetization M p 
in the presence of both Zccman and orbital effects is given 
by^ 



M p = M n [l + f(H)} = M n + M de 



(2) 



where f(H) is a field dependent function and Mdev is the 
deviation of the mixed-state paramagnetic magnetization 
from the linear in H behavior, i.e. from M n . Therefore, 
in order to determine M p , hence Mdev, one needs to de- 
termine f(H). We determine f(H) from torque measure- 
ments in the mixed and normal states, as follows. 

The magnetic moment of a sample placed in a mag- 
netic field feels a torque f = M x H . Hence, both the 
paramagnetic and vortex magnetizations have associated 
induced torques r p and t v , respectively. As we have pre- 
viously shown^ the reversible torque measured in the 
mixed state is given by: 



T rev (T, H, 9) =T p + T v 



(3) 



where 



Tp (T,H,0) = T n [l + f(H)] 



Xl{T)- X c n {T) 



H 2 sin 20[1 + f{H)\ = A(T, H) sin 29, 



(4) 

with A(T, H) a fitting parameter, and t v is given by Ko- 
gan's models Equation (4) is valid if the magnetizations 
and M£ along the a and c crystallographic direc- 
tions, respectively, have the same H dependence; i.e., if 
the function f(H) is direction independent. f(H) can 
then be obtained from Eq. (4) as 



in which A(T, H) and B 



(5) 



H 2 are obtained by 



fitting the torque data in the mixed and normal state, 
respectively. Note that f(H) — in the normal state due 
to the fact that M p = M n . Therefore, knowing f(H), one 
can obtain the mixed state paramagnetic magnetization 
M p from Eq. (2) and the vortex magnetization M v by 
subtracting M p from the measured magnetization in the 
mixed state [see Eq. (1)]. 

We note that the above assumption that the depen- 
dence of f(H) on direction is negligible is supported by 
the present torque data, which can be fitted only with 
a A(T, H) sin 29 dependence, with no additional angular 
dependences. This assumption that the magnetizations 
Mp and Mp along a and c axis have the same H depen- 
dence is, in addition, supported by previous studies. For 



example, we have previously shown 2 ^ that the field de- 
pendent in-plane normal-state resistivity data measured 
along the c and a crystallographic directions scale, with 
the anisotropy as the scaling factor (see Fig. 3 of the 
above reference) . This implies that the same field depen- 
dence, hence same physics, dominates the charge trans- 
port when H is applied along the a and c directions. In 
another study, 27 which points toward the same conclu- 
sion, the authors have shown that the difference between 
the response of a high temperature superconductor in the 
mixed state when the magnetic field is along the a and c 
directions is closely related with the field dependence of 
the upper critical field. In fact, the authors have shown, 
through calculations of thermodynamic and electromag- 
netic properties, the presence of a similar scaling law for 
several thermodynamic properties, including magnetiza- 
tion. Since the high temperature superconductors are 
generally even more anisotropic than CeCoIns , we believe 
that these results most likely apply also to this latter sys- 
tem; i.e. the spin scattering along the c and a directions 
of CeCoIn 5 has the same field dependence, but different 
coefficients, which are related with the anisotropy. 

We performed torque measurements on CeCoIn 5 single 
crystals both in the normal and mixed states. From nor- 
mal state torque measurements we obtain B — (2.39 x 
10 _7 if 2 ) Nm, where H is in Tesla. We have already 
show»i£ that t v and t p can be successfully separated in 
the mixed state, with t v well described by Kogan's model 
and T p — .A sin 20 [see Eqs. (3) and (4)]. Therefore, we 
obtain A(H), shown in the inset to Fig. 3(a), by fit- 
ting the torque data in the mixed state with Eq. (3). 
A simple fit of these A(H) data with a power law gives 
A(H) = (1.57 x 10 _7 .ff 2 - 32 ) Nm, where H is in Tesla 
(the solid line in the inset). Note that the magnetic field 
dependence of the coefficient A, which gives the field de- 
pendence of the paramagnetic contribution in the mixed 
state, is stronger than H 2 , which is typical for param- 
agnetism. Also, we note that the plot of A(H) has data 
points only up to 1.8 T since Kogan's model, which gives 
the vortex torque, is valid only for fields much smaller 
than (see Ref. 18 for more discussion). 

The plot of f(H) for the field range 0-18 kOe, ob- 
tained from Eq. (5), is shown in the main panel of Fig. 
3(a). As discussed above, knowing f(H), one can ob- 
tain the paramagnetic and vortex magnetizations in the 
mixed state. Figure 3(b) shows the field dependence of 
different magnetization curves. The diamonds give the 
vortex response M v , obtained by subtracting M p [given 
by Eq. (2)] from Eq. (1). Since the analysis of the torque 
data is limited to magnetic fields lower than ~ 18 kOe^ 
there are no data points in M V (H) in the field region close 
to ■ However, a linear extrapolation of the available 

high field data leads exactly to H]}£ [see the dashed line 
in Fig. 3(b)]. This linear extrapolation of M V (H) is rea- 
sonable since the vortex magnetization should be linear 



in H when H is close to H. 



c2 ' 



Knowing M V (H) up to Hj^, permits the calculation of 
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f(H) for H > 18 kOe [see the data points for H > 20 
kOeinFig. 3(a)] from Eq. (2) with M p given by Eq. (1). 
Finally, knowing f(H) over the whole H range allows the 
determination, from Eq. (2), of Mdev(H), shown by the 
reversed solid triangles in Fig. 3(b), and M p (H) shown 
in Fig. 4. Note that Mdev is a non-monotonic function 
of H. The shapes of Md ev {H) and Mi(H) (open circles) 
are similar, which shows that the anomalous behavior of 
Mi{H) is due to Md ev {H). This reinforces the suitability 
of our analysis. Also note that M p (H) is not linear in H 
in the mixed state. 

The inset to Fig. 4 is a plot of the differential para- 
magnetic susceptibility Xp = dM p /dH. Note that the 
paramagnetic susceptibility in the mixed state is mag- 
netic field dependent, while its value is constant, equal 
with 1.84 x 10 -5 emu/g in the normal state. The jump 

in Xp(H) around reflects the superconducting phase 
transition. 

The M p {H) and XpC^O dependences below H]}£ can be 
understood from the fact that the "heavy" electrons of 
CeCoIns contribute to both paramagnetism and super- 
conductivity, and the Zeeman effect is large. Specifically, 
in the normal state, the large paramagnetic moment 
comes from the heavy fermion quasi-particles. In the 
mixed state, the condensation energy favors the forma- 
tion of Cooper pairs with one spin up and one spin down, 
(note that CeCoIns has a <i-wave symmetry, i.e., singlet 
spin pairing), while the large Zeeman effect decouples 
the spins of some of the "heavy" electron Cooper pairs, 
which, hence, contribute to paramagnetism. Therefore, 
the magnetization M p (H) in the mixed state is sup- 
pressed compared with the magnetization M n (H) in the 
normal state. The number of decoupled " heavy" electron 
spins available to align with H, hence to participate in 
the mixed state paramagnetism, increases with increas- 
ing H. This gives rise to an increase in \ p with increasing 
H. The finite value of Xp(H) as H — > is consistent with 
the finite density of quasiparticles present in the mixed 
state. As expected, this value of \p is smaller than the 
value in the normal state, and it is very close to the value 
reported by NMR measurements £ 

In materials in which the electrons responsible for para- 
magnetism do not participate in superconductivity (e.g., 
localized d or / electrons), the susceptibility in the mixed 
state is field independent and hence the paramagnetic 
magnetization is a linear extrapolation of the normal 
state paramagnetism. Our result of a suppressed para- 
magnetism in the mixed state is consistent with recent 



115 In and 59 Co NMR measurements £ 

Note that the vortex response in the mixed state has 
a monotonic field dependence [see M V (H) in Fig. 3(b)] 
with no anomaly observed for the investigated tempera- 
ture T = 1.8 K {T/T c = 0.78). Theorectical calculations 
of Adachi et alj^ have shown an anomalous response, 
i.e. a change in the M v vs H curvature below the re- 
duced temperature T/T c = 0.3 with no anomaly above 
this reduced temperature. Hence, our experimental re- 
sult confirms this latter theoretical prediction. 

We should mention that we also measured M(H) for 
H || a-axis. However, in the temperature and field range 
investigated (T > 1.76 K and H < 50 kOe), no second 
peak was obtained after subtracting the linear param- 
agnetic moment (see upper inset to Fig. 2). Neverthe- 
less, note that this M\(H) curve is still anomalous in 
the sense that there is a change of curvature, which im- 
plies that a similar anomalous paramagnetism exists in 
the H || a direction due to the presence of large Zee- 
man effect. However, the larger upper critical field along 
the a-axis requires even lower temperatures and higher 
magnetic fields for the full observation of this effect. 

D. Summary 

In summary, we performed magnetization and torque 
measurements both in the normal and mixed states of 
CeCoIns single crystals in order to study the param- 
agnetic and vortex response in the presence of a large 
Zeeman effect present in this material. The paramag- 
netic magnetization is suppressed in the mixed state 
and the spin susceptibility is field dependent, increasing 
with increasing field. This H dependence is a result 
of the fact that heavy electrons contribute to both 
superconductivity and paramagnetism and the Zeeman 
effect is large in this material. There is no anomaly 
present in the vortex response in the temperature range 
investigated. 
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I. FIGURE CAPTIONS 

Figure 1. (Color online) Magnetic field H dependence 
of the dc magnetization M mes measured at 1.76 K with 
H || c-axis on a CeCoIns single crystal. The solid line is 
a linear fit of M mes (H) in the normal state. Inset: Upper 

lie 

critical field parallel to the c-axis HJ 2 — temperature T 
phase diagram. The open squares are data taken from 
RefJ£ while open circles are data extracted from present 
M mes (H) measurements. 

Figure 2. (Color online) Magnetic field H dependence 
of the magnetization Mi measured at 1.76 K which is 
obtained by subtracting the paramagnetic contribution 
as an extrapolation of the normal state paramagnetism. 
Lower inset: Plot of Mi(H) measured at 1.76, 1.80, 1.85, 
1.90, 1.95, 2.00, 2.05, and 2.10 K. Upper inset: Mag- 
netic field H dependence of the dc magnetization M mes 
measured at 2 K for H || a. 

Figure 3. (Color online) (a) Plot of field H depen- 
dence of the function / determined at 1.8 K. The solid 
line is a guide to the eye. Inset: H dependence of the 
fitting parameter A. The solid line is a fit of the data 
with a simple power law. (b) H dependence of vortex 
magnetization M v (solid diamonds), deviation magneti- 
zation Mdev (solid reversed triangles) , and magnetization 
Mi data of Fig. 2 (open circles) of CeCoIns measured at 
1.8 K. The dashed line in M V (H) is a linear extrapola- 
tion of the high field data. The solid lines in M v (H) and 
Mdev{H) are guides to the eye. 

Figure 4. Magnetic field H dependence of the para- 
magnetic magnetization M p . Inset: H dependence of 
differential susceptibility x = dM/dH. The solid line is 
a guide to the eye. 
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